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(57) Abstract: A Doppler spread for a communications channel is 
measured by providing an estimate of the com muni cations channel 
and generating an autocorrelation function for the estimate of the 
communications channel. One of a plurality of autocorrelation function 
hypotheses is selected to approximate the autocorrelation function for the 
estimate of the communications channel wherein each of the autocorrelation 
function hypotheses corresponds to a respective Doppler spread estimate 
hypothesis. One of the Doppler spread estimate hypotheses is selected 
corresponding to the selected autocorrelation function hypotheses as an 
estimate of the Doppler spread for the communications channel. 
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DOPPLER SPREAD ESTIMATION USING CHANNEL AUTOCORRELATION FUNCTION HYPOTHESES 



FIELD OF THE INVENTION 
The present invention relates to the field of communications and more 
particularly to receiving radio communications. 

5 BACKGROUND OF THE INVENTION 

A radio channel for a mobile terminal in a cellular radiotelephone 
communications system may be difficult to operate. In particular, the 
transmitted signals are often reflected, scattered, diffracted, delayed, and 
attenuated by the surrounding environment. Moreover, the radio channel for 

10 a mobile terminal is often not stationary because of movement of the mobile 
terminal and movement of objects near the mobile terminal. The mobile 
terminal may move rapidly when used in an automobile, and other vehicles 
may also be in motion near the mobile terminal. 

Characteristics of the radio channel may also vary from one area to 

1 5 another due to differences in terrain/buildings, climate, and/or other factors. 
The propagation of a radio signal along the radio channel may thus be subject 
to multi-path fading, shadowing, and path loss. Of these factors, multi-path 
fading may be the most significant, and multi-path fading can be characterized 
by envelope fading, Doppler spread, and time-delay spread. 

20 Doppler shift is the frequency shift experienced by the radio signal 

when the mobile terminal is in motion, and the Doppler spread is a measure of 
the spectral broadening caused by the time rate of change of the mobile radio 
channel. Doppler spread thus leads to frequency dispersion, and the Doppler 
spread in the frequency domain is closely related to the rate of change in the 

25 observed signal. The adaptation time of an algorithm used in an adaptive 
receiver should thus be faster than the rate of change of the channel to be 
able to accurately track the fluctuations in the received signal. 

A mobile terminal in a DAMPS cellular radiotelephone communications 
system, for example, may experience a Doppler spread in the range of 0Hz to 

1 
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250Hz depending on the vehicle speed, carrier frequency, and other factors. 
Knowledge of the rate of change of the radio channel can be used to improve 
receiver performance and/or reduce receiver complexity. Moreover, the 
adaptation parameters for an adaptive receiver can be varied as a function of 
5 the Doppler spread. Instead of fixing the tracker and interpolation parameters 
for the worst case expected Doppler spread, for example, the parameters can 
be changed adaptively as a function of the Doppler information to improve 
performance. Similarly, Doppler spread information can be used to control 
the receiver adaptively for different speeds at which the mobile may be 

10 traveling. In other words, different receiver algorithms can be used depending 
on the speed at which the mobile terminal is traveling. 

Estimates of the Doppler spread can thus be used to increase receiver 
performance. The parameters of a receiver adaptation algorithm can be 
varied as a function of the Doppler spread to adaptively optimize a coherent 

15 detector in a receiver, for example. In addition, the hand-off process in a 
cellular mobile telephone system can be enhanced if an estimate of the 
Doppler spread is available. Handoff of a fast moving mobile terminal to a 
micro cell can thus be avoided. 

Doppler spread estimation is discussed, for example, in U.S. Patent 

20 No. 4,723,303 to Koch entitled "METHOD OF AND CIRCUIT 

ARRANGEMENT FOR MEASURING THE QUALITY OF RADIO- 
TRANSMISSION CHANNELS OF A RADIO-TRANSMISSION SYSTEM", and 
U.S. Patent No. 5,016,017 to Raith entitled "METHOD OF CONTROLLING 
THE FREQUENCY OF A COHERENT RADIO RECEIVER AND APPARATUS 

25 FOR CARRYING OUT THE METHOD". The disclosures of each of these 
patents is hereby incorporated herein in their entirety by reference. 

A method of estimating Doppler spreads from a sequence of channel 
estimates is discussed, for example, in the thesis by Lars Lindbom entitled 
"Adaptive Equalization For Fading Mobile Radio Channels", (Techn.Lic.Thesis 

30 No. UPTEC 921 24R, November 1992, Department of Technology, Uppsala 
University, Uppsala Sweden), the disclosure of which is hereby incorporated 
herein in its entirety by reference. In the Lindbom thesis, differentials of the 
channel estimate, which comprise differences of values between two points in 
time, are used to estimate Doppler spreads. These differentials, however, 

2 
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may be noisy so that averaging may be needed. As a result, the averaging 
may give biased estimates of the Doppler spread. 

Other Doppler spread estimation techniques are discussed in the 
reference by Karim Jamal et al. entitled "Adaptive MSLE Performance On The 
5 D-AMPS 1900 Channel" (IEEE Trans. Vehic. Technol., vol. 46, Aug. 1997), 
and the reference by M. Morelli et al. entitled "Further Results In Carrier 
Frequency Estimation For Transmissions Over Flat Fading Channels" (IEEE 
Commun. Letters, vol. 2, pp. 327-330, Dec. 1998). The disclosures of these 
references are also incorporated herein in their entirety by reference. 
10 Notwithstanding the approaches discussed above, there continues to 

exist a need in the art for improved Doppler spread estimation approaches. 

SUMMARY OF THE INVENTION 
It is therefore an object of the present invention to provide improved 
1 5 methods of estimating Doppler spreads for communications channels and 
related systems and receivers. 

It is another object of the present invention to provide less complex 
methods of estimating Doppler spreads and related systems and receivers. 
These and other objects can be provided according to the present 
20 invention by providing an estimate of the communications channel and 

generating an autocorrelation function for the estimate of the communications 
channel. One of a plurality of autocorrelation function hypotheses is selected 
to approximate the autocorrelation function for the estimate of the 
communications channel wherein each of the autocorrelation function 
25 hypotheses corresponds to a respective Doppler spread estimate hypothesis. 
The Doppler spread estimate hypothesis corresponding to the selected 
autocorrelation function hypothesis is then selected as an estimate of the 
Doppler spread for the communications channel. 

The autocorrelation function hypotheses can thus be saved in a 
30 memory of a Doppler spread estimator according to the present invention and 
compared to the autocorrelation function of the estimate of the 
communications channel, with the closest autocorrelation function hypotheses 
being accepted as an estimate of the actual autocorrelation function for the 
communications channel. The Doppler spread hypothesis corresponding to 

3 
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the closest autocorrelation function can thus be used as an estimate of the 
actual Doppler spread for the communications channel. The complexity of 
calculations used to estimate the Doppler spread can thus be reduced while 
providing a relatively accurate estimate of the Doppler spread. In addition, the 
5 number of autocorrelation function hypotheses used can be increased to 
provide a more accurate estimation, or decreased to reduce the number of 
calculations and amount of memory used. 

In greater detail, the selection of one of the autoconnection function 
hypotheses can include comparing the autocorrelation function for the 

10 estimate of the communications channel with each of the plurality of 

autocorrelation function hypotheses. In addition, the selection of one of the 
autocorrelation function hypotheses can include selecting one of the plurality 
of autocorrelation function hypotheses most closely approximating the 
autocorrelation function for the estimate of the communications channel. 

1 5 Moreover, the selection of one of the autocorrelation function 

hypotheses can include generating a plurality of error signals respectively 
corresponding to the plurality of autocorrelation function hypotheses wherein 
each of the error signals represents a difference between the respective 
autocorrelation function hypothesis and the autocorrelation function of the 

20 estimate of the communications channel, and comparing the error signals to 
select the autocorrelation hypothesis to approximate the autocorrelation 
function for the estimate of the communications channel. In particular, the 
error signals can be compared to choose the error signal representing a least 
difference between the corresponding autocorrelation function hypothesis and 

25 the autocorrelation function for the estimate of the communications channel. 
In addition, the comparison of the error signals can be preceded by averaging 
each of the error signals to provide averaged error signals wherein comparing 
the error signals comprises comparing the averaged error signals. 

Methods, systems, and receivers according to the present invention 

30 can thus be used to provide estimates of Doppler spreads for a 
communications channel with reduced complexity. 



4 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a block diagram of a communications system including a 
receiver according to the present invention. 

Figures 2-4 are block diagrams of receivers according to the present 

5 invention- 
Figure 5 is a block diagram of a Doppler spread estimator according to 

the present invention. 

Figure 6 is a graph illustrating autocorrelation functions for radio 
channels at different speeds of a receiver relative to a base station. 
10 Figure 7 is a table illustrating a storage of samples of the 

autocorrelation functions of Figure 6. 

DETAILED DESCRIPTION 
The present invention will now be described more fully hereinafter with 
15 reference to the accompanying drawings, in which preferred embodiments of 
the invention are shown. This invention may, however, be embodied in many 
different forms and should not be construed as limited to the embodiments set 
forth herein; rather, these embodiments are provided so that this disclosure 
will be thorough and complete, and will fully convey the scope of the invention 
20 to those skilled in the art. Like numbers refer to like elements throughout. 

Figure 1 illustrates a transmitter T and a receiver R according to the 
present invention wherein data d is transmitted by the transmitter T over the 
radio channel c. The received signal r is a function of the transmitted data d, 
the radio channel c. and noise n. In a flat fading channel: 

25 

r=h»d+n. (equation 1) 

As discussed above, receiver performance can be improved by estimating the 
Doppler spread and using the estimated Doppler spread to adapt receiver 
30 functions. More particularty, the estimated Doppler spread can be used to 
more accurately estimate the radio channel h. The use of Doppler spread 
estimators in receivers is discussed in co-pending U.S. Patent Application 
Serial No. .to Leonid Krasny et al. entitled "DOPPLER SPREAD 
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ESTIMATION SYSTEM" filed on March 22, 1999 (Docket No. P10367-RCUR 
and 1280.00105). The Krasny et al. application is assigned to the assignee of 
the present application, and the Krasny et al. application shares a common 
inventor with the present application, in addition, the disclosure of the Krasny 
5 et al. application is hereby incorporated herein in its entirety by reference. 

Various receivers Ra, Rb, and Rc including Doppler spread estimators 
according to the present invention are illustrated in Figures 2-4. In particular, 
the receiver Ra of Figure 2 is adapted for use with known pilot symbols, and 
this receiver includes an antenna 21 that receives the signal r, a radio 

10 receiver and converter 23a, a channel estimator 25a, a Doppler spread 

estimator 27a, a known symbol block 29a, and a signal processor 31a. The 
antenna 21a receives the radio signals r, and the radio receiver and converter 
23a filters, amplifies, and converts the received radio signals r into digital 
samples for processing. More particularly, the received radio signals can be 

15 converted to a form such as complex numerical sample values suitable for 
processing. The channel estimator receives the converted radio signals and 
known symbols from the known symbol block 29a. 

In particular, the known symbol block 29a can provide pilot symbols or 
other known reference symbols, such as synchronization symbols, included in 

20 the received signal r that can be used to calculate channel estimates. 

Moreover, the known symbol block 29a can include a memory where the 
known symbols are stored or a code generator that can generate the known 
symbols. The channel estimator 25a correlates the received digital samples 
with the known symbols to provide an estimate of the channel c to be used by 

25 the Doppler spread estimator 27a. The Doppler spread estimator estimates 
the Doppler spread using the channel estimates and then sends the estimated 
Doppler spread to the signal processor 31a. The operation of the Doppler 
spread estimator will be discussed in greater detail below. The signal 
processor 31a processes the sampled signal to extract information, and the 

30 signal processor 31a provides feedback to the channel estimator 25a so that 
the channel estimates can be improved after the Doppler spread estimation. 

The receiver Rb of Figure 3 is adapted for use without known pilot 
symbols. This receiver includes an antenna 21b that receives the signal r, a 
radio receiver and converter 23b, a channel estimator 25b, a Doppler spread 



BNSDOCID: <WO 0113537A1 I > 



WO 01/13537 PCT/US00/21081 

estimator 27b, a symbol estimator 29b t and a signal processor 31b. The 
receiver Rb is similar to the receiver Ra of Figure 2 with the exception that the 
symbol estimator 29b is used instead of the known symbol block 29a used in 
the receiver of Figure 2. The symbol estimator 29b can be used in 
5 applications where symbols are not known. The symbol estimator 29b, for 
example, can use error detection and correction techniques to estimate 
received symbols with a high level of accuracy. These estimated symbols can 
then be used by the channel estimator 25b to estimate the channel h. The 
Doppler spread estimator 27b uses the channel estimates to estimate the 

10 Doppler spread as discussed in greater detail below. 

In code-division multiple access (CDMA) cellular systems (such as IS- 
95 systems), a transmitter transmits a stream of known symbols known as the 
pilot code. The pilot code is transmitted on the same channel and at the 
same time as other information bearing symbols using different spreading 

15 codes. The receiver Rc of Figure 4 provides Doppler spread estimations in 
such a CDMA system. The receiver Rc of Figure 4 is adapted for use with 
known pilot symbols, and this receiver includes an antenna 21c that receives 
the signal r, a radio receiver and converter 23c, a channel estimator 25c, a 
Doppler spread estimator 27c, and a signal processor 31c. In this CDMA 

20 receiver, the channel can be estimated directly by the channel estimator 25c 
without the known symbol block 29a of Figure 1 or the symbol estimator 29b 
of Figure 2, and the channel estimates used for Doppler spread estimation. 

The channel estimator 25c correlates the received signal r including 
the known pilot code and other codes in additive superimposition, and filters 

25 the resultant complex correlation to obtain channel estimates. The received 
signal can also be correlated with other codes carrying information to be 
decoded. The results of correlating with information carrying codes are 
multiplied by the conjugate of pilot code correlations for the same delay, and 
the results added to coherently combine the multi-path signals. In wide band 

30 CDMA (WBCDMA) systems, modulation symbol intervals may be much 
shorter thus allowing multiple propagation paths to be resolved with much 
finer time resolution. 

The receivers of Figures 2-4 thus illustrate various receivers including 
Doppler spread estimators according to the present invention. In each 



BNSDOCID: <WO 0113537A1_U> 



WO 01/13537 PCT/USOO/21081 

receiver, channel estimates are provided to the Doppler spread estimator for 
calculation of the Doppler spread estimates. Receivers including Doppler 
spread estimators according to the present invention are not limited to the 
channel estimators discussed above, and those having skill in the art will 
5 recognize that any channel estimation technique can be used to provide 
channel estimates. 

In particular, the channel estimator can generate an estimate of the 
channel over a TDMA time slot using symbols representing data samples 
received during the time slot. For example, the receiver Rb of Figure 3 can 

10 include the symbol estimator 296 to estimate the symbols. The Doppler 

spread estimator can then use the channel estimate for the current time slot to 
calculate a Doppler spread estimate for the current time slot, and the signal 
processor can use this Doppler spread estimate in calculations for 
subsequent time slots, such as calculations of channel estimates for 

1 5 subsequent time slots. In other words, the Doppler spread estimate for the 
current time slot can be used to update a long term Doppler estimate in the 
signal processor which is used in subsequent calculations. The long term 
Doppler estimate can be updated, for example, using averaging techniques 
discussed in greater detail below. 

20 When a symbol estimator is used as discussed with reference to Figure 

3, it may be useful to perform a cyclical redundancy check (CRC) on the 
estimated symbols used to calculate the Doppler spread estimate for the time 
slot. If the cyclical redundancy check passes, the channel estimate should be 
relatively accurate, and the resulting Doppler spread estimate can be used to 

25 update the long term Doppler estimate. If the cyclical redundancy check fails, 
however, the channel estimate may be unreliable, and it may thus be 
desirable to not update the long term Doppler estimate using a Doppler 
spread estimate based on a potentially unreliable channel estimate. 

The Doppler spread estimation will now be discussed in greater detail 

30 with reference to the Doppler spread estimator of Figure 5. The Doppler 

spread estimator can be used, for example, with receivers of radiotelephone 
communication systems operating according to either the DAMPS or 
DAMPS+ standards. In DAMPS+, known pilot symbols are provided. In 
DAMPS, with a multi-pass demodulation, the Class 1 bits demodulated in the 
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first pass can be used as pilot symbols if they pass the cyclical redundancy 
check (CRC). Preferably, phase ambiguities between the pilot slots is 
reduced to increase the reliability of the results. A method of resolving phase 
ambiguity is discussed, for example, in the reference by T. Fulghum entitled 
5 "Channel Interpolation On Second Pass Demodulation", (Tech. Rep. Tr/X 
98:1230, Ericsson, RTP, NC, Feb. 22, 1999) the disclosure of which is hereby 
incorporated herein in its entirety by reference. 

Assuming that the pilot symbols are available and that the phase 
ambiguity between pilot symbols has been reduced to an acceptable level, 

10 samples of the autocorrelation function R e can be found by the 

autocorrelation calculator 51 using the channel estimates obtained from the 
pilot symbols. The samples of the autocorrelation function are then compared 
with different hypotheses (Hi, H 2 , H 3 , ... H k ) of the true autocorrelation 
function of different Doppler spread values using the hypothesis comparator 

15 53 to provide error calculations (e 1(m) , e 2 (m), e3(m), -» Ckm) for the 

autocorrelation function samples with respect to each of the hypotheses. 
Each of the error calculations can then be averaged over several slots using 
the averager 55 to provide a respective plurality of averaged error values 
(e a y.f, em e av , 3 , ... e ay , k ). The Doppler spread corresponding to the 

20 hypothesis resulting in the lowest averaged error value is then selected by the 
minimum error selector 57 to provide the estimated Doppler spread fdoppier- 

In greater detail, the channel estimates can be determined using 
techniques available now or in the future, and the hypotheses of the different 
autocorrelation functions (H^ H 2 , H 3 , ... H k ) corresponding to different Doppler 

25 spread values can be calculated and stored in memory. In particular, actual 
radio channels and corresponding Doppler spread values can be determined 
at different mobile terminal speeds relative to the base station, and the 
resulting autocorrelation functions can be calculated for each speed. 
Graphical examples of hypotheses of different autocorrelation functions 

30 (correlation vs. x) corresponding to radio channels measured at different 

speeds and corresponding to different Doppler spread values are illustrated in 
Figure 6. 
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Samples of each autocorrelation function representing a respective 
hypothesis Hr can be stored in memory for the hypothesis comparator 53 to 
provide the hypothesis Hi, H 2( H 3 , ... Hk shown in Figure 5. For example, 
samples of each autocorrelation function of the graph of Figure 6 
5 corresponding to different values of t can be stored in memory as illustrated 
in Figure 7. While five different hypotheses are shown in the graph of Figure 
6, any number can be used. A greater number (k) of hypotheses can be used 
to provide a greater accuracy in the estimation of Doppler spread values, 
while a fewer number (k) of hypotheses can provide less complicated 

10 operation with fewer calculations and a lower memory requirement. 

Moreover, the number of samples z saved for each hypothesis can be varied 
depending on the desired level of accuracy for the Doppler spread estimation. 

The autocorrelation function estimation for a channel with frequency 
error can be calculated using the autocorrelation estimator 51 using the 

15 channel estimates over the known field as: 

* a (c) = -^— Y J k(k)hl(k~x). (equation2) 

In this equation, R e is the estimated autocorrelation function including the 
20 frequency error, and f? 0 is the channel estimates over the known fields (or pilot 
symbols). Note that when there is a frequency error in the received signal, it 
may directly effect the channel estimates which may be rotated by the amount 
of the frequency error thus effecting the autocorrelation estimate. Designating 
the frequency error as f e , the rotation in the channel estimates can be 
25 calculated as: 

h e (k) = he jWT , (equations) 

where h represents the channel estimates when there is no frequency error. 
30 The autocorrelation function with frequency error can thus be rewritten as: 
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15 



25 



R.b) = -T- f J kk)e J2 ^ T ^(k-'z)^ k -^, (eq uation4) 

L—T 



or: 



R e (T) = R(X )e i2Tf * T , (equation 5) 

5 

where R (t) is the autocorrelation estimate without the frequency error. 

The frequency error f e can be calculated using know methods as 
discussed for example in the reference by Morelli et al. entitled "Further 
Results In Carrier Frequency Estimation For Transmissions Over Flat Fading 
10 Channels", (IEEE Commun. Letters, vol. 2. pp. 327-330, Dec. 1998), the 

disclosure of which is hereby incorporated herein in its entirety by reference. 
The effect of the frequency error can thus be removed from the correlation of 
the channel estimates as: 



R(x ) — R e ( T )e J2 *' zT , (equation 6) 



The frequency error-free correlation estimates can then be compared with the 
different hypotheses (Hi, H 2 , H 3 . ... H k ) of the hypothesis comparator 53. 
The hypotheses (Hi, H 2 , H 3 , ... H K ) can be determined and stored in 
20 memory as discussed above with regard to Figures 6 and 7. The error e k (m) 
of the autocorrelation estimate R e with respect to each of the hypotheses is 
calculated using: 



e t {m) = ^R(n)-R(P)R HJl {nf. (equation 7) 

#1=1 

In this equation, e k (m) is the error corresponding to the mth slot between the 
estimated autocorrelation function and the different samples of the /cth 
hypothesis of the normalized true autocorrelation function RhAh) given as: 
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(equation 8) 



where fdoppter.k is the /cth Doppier spread hypothesis. The error term e k is 
averaged over several slots at the averager 55 (filter) using the equation: 

e <»* = TtX**^' (equation 9) 



where e QVtk is the averaged error corresponding to the /cth hypothesis, and M 

is the length of the averaging window. In this example, block averaging is 
10 used to average the error. Alternately, other averaging techniques such as 

running averaging or sliding window averaging can be used! Averaging thus 

reduces statistical errors. 

The minimum error selector 57 determines the lowest of the averaged 

errors e aVfk to select the autocorrelation function hypothesis H k most closely 
15 approximating the autocorrelation function for the received channel. The 

Doppler spread hypothesis selector selects the Doppler spread hypothesis 

corresponding to the selected autocorrelation function hypothesis as the 

estimate of the Doppler spread / do PP /er- In other words, the minimum of the 
averaged errors e aVt i-e aVtk is used to select the autocorrelation function 
20 hypothesis resulting in the lowest averaged error (at the minimum error 
selector 57), and the Doppler spread hypothesis corresponding to the 
selected autocorrelation function hypothesis is selected as the estimate of the 
Doppler spread (at the Doppler spread hypothesis selector 59) according to 
the following equation. 

25 

■A 

doppler & wmn ^ov,** (equation 10) 

J doppier,k 

Alternately, the autocorrelation function with frequency error can be 
calculated as discussed in the Morelli et al. reference. Instead of finding the 
30 frequency error and removing it from the estimated autocorrelation function, 
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we can take the absolute value of the estimated autocorrelation to obtain the 
envelope autocorrelation function as: 



R(t) = \R (x)\ = | A {x)e IZ7lf e TT 1 = 1* (equation 1 1 ) 

5 

This absolute value of the autocorrelation function estimate (without frequency 
error) can then be compared with the different hypotheses of the true 
envelope autocorrelation function corresponding to the different Doppler 
spread values as: 

10 

e k (m) = YA R ( n )- R (°) R H*( n A (equation 12) 

where, 

15 = \ J o (Appier* nT)\ . (equation 1 3) 

The error terms are used as before to determine the hypothesis that provides 
the least error in the estimation of the Doppler spread value. 

In DAMPS systems, phase ambiguity can be substantially removed 

20 within a time slot using known techniques. Phase ambiguity across time slots, 
however, may still exist. According to the present invention, the correlation 
within slots can be estimated and compared with the different hypotheses. 
Errors across slots, however, are averaged. 

As previously discussed, the number of hypotheses can be varied 

25 depending on the desired accuracy and complexity of the Doppler spread 

estimator and receiver. A greater number of hypotheses can provide greater 
accuracy while requiring more memory and computation. A fewer number of 
hypotheses can reduce memory use and computation while providing a less 
accurate Doppler spread estimation. 

30 The Doppler spread estimation methods according to the present 

invention can thus be implemented with relative computational efficiency, and 
good results. In particular, the use of the envelope correlation function can 
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provide a level of performance close to that of other techniques without 
requiring frequency error estimation. 

Doppler estimation, according to the present invention, can be 
performed, for example, in the context of a DAMPS cellular radio telephone 
5 system using the following downlink slot and transmission formats: tt/4- 
DQPSK with root raised cosine pulse shaping; TDMA frame of 20ms; and 3 
users sharing the TDMA frame. Moreover, each user can transmit twice 
during a frame with a slot duration of 6.667 ms. The transmission medium 
can be a Rayleigh fading channel, and the channel can be simulated using 

10 Jake's fading model. A single antenna receiver and a carrier frequency of 
1900 MHz can be used. 

In downlink joint demodulation, only one Doppler spread needs to be 
estimated for both desired and interfering signals because both desired and 
interfering base stations are fixed and the mobile terminal is moving. This 

15 assumption is reasonable assuming there is no direct path. The Doppler 
spreads can thus be estimated using information corresponding to the 
stronger user which is the desired user in most cases. 

In the receiver, ideal sync location can be assumed to be known, and 
the received samples can be used to coherently demodulate the user's 

20 information sequences. Kalman based channel estimation using a second 
order autoregressive model (AR-2) can be used for channel tracking. As 
discussed above, channel tracking depends on the Doppler information. As 
the Doppler information may not be known initially, the tracker parameters can 
be initially set to a relatively high Doppler spread of 100Hz. As the Doppler 

25 spread estimate converges to the actual Doppler spread, the tracker 
parameters should also be changed in the adaptive receiver. 

The present invention may be embodied as methods or devices. In 
addition, the invention may take the form of an entirely hardware embodiment, 
an entirely software embodiment, or an embodiment combining both hardware 

30 and software aspects. The present invention has been described in part with 
respect to the block diagrams of Figures 1-5. It will be understood that each 
block of the illustrations, and combinations of blocks, can be implemented by 
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computer program instructions. These program instructions, which may 
represent steps, may be provided to a processor to produce a machine. 

Accordingly, blocks of the block diagrams support combinations of 
means for performing the specified functions in combinations of steps for 
5 performing the specified functions. It will be understood that each block of the 
illustrations, and combinations of blocks, can be implemented by special 
purpose hardware-based systems which perform the specified functions or 
steps, or combinations of special purpose hardware and computer 
instructions. 

10 In the drawings and specification, there have been disclosed typical 

preferred embodiments of the invention and, although specific terms are 
employed, they are used in a generic and descriptive sense only and not for 
purposes of limitation, the scope of the invention being set forth in the 
following claims. 
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THAT WHICH IS CLAIMED IS: 

1 . A method for estimating a Doppler spread for a communications 
channel, the method comprising: 

providing an estimate of the communications channel; 
generating an autocorrelation function for the estimate of the 
5 communications channel; 

selecting one of a plurality of autocorrelation function hypotheses to 
approximate the autocorrelation function for the estimate of the 
communications channel wherein each of the autocorrelation function 
hypotheses corresponds to a respective Doppler spread estimate hypothesis; 
10 and 

selecting one of the Doppler spread estimate hypotheses 
corresponding to the selected autocorrelation function hypotheses as an 
estimate of the Doppler spread for the communications channel. 

2. A method according to Claim 1 wherein selecting one of the 
autocorrelection function hypotheses comprises comparing the 
autocorrelation function for the estimate of the communications channel with 
each of the plurality of autocorrelation function hypotheses. 

3. A method according to Claim 1 wherein selecting one of the 
autocorrelation function hypotheses comprises selecting one of the plurality of 
autocorrelation function hypotheses most closely approximating the 
autocorrelation function for the estimate of the communications channel. 

4. A method according to Claim 1 wherein selecting one of the 
autocorrelation function hypotheses comprises: 

generating a plurality of error signals respectively corresponding to the 
plurality of autocorrelation function hypotheses wherein each of the error 
5 signals represents a difference between the respective autocorrelation 
function hypothesis and the autocorrelation function of the estimate of the 
communications channel; and 
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comparing the error signals to select the autocorrelation hypothesis to 
approximate the autocorrelation function for the estimate of the 
10 communications channel. 

5. A method according to Claim 4 wherein comparing the error 
signals comprises choosing an error signal representing a least difference 
between the corresponding autocorrelation function hypothesis and the 
autocorrelation function for the estimate of the communications channel. 

6. A method according to Claim 5 wherein choosing the error 
signal comprises choosing a least one of the error signals. 

7. A method according to Claim 4 wherein comparing the error 
signals is preceded by: 

averaging each of the error signals to provide averaged error signals 
wherein comparing the error signals comprises comparing the averaged error 
5 signals. 

8. A method according to Claim 7 wherein averaging each of the 
error signals comprises one of block averaging, running averaging, and sliding 
window averaging. 

9. A method according to Claim 1 wherein each of the plurality of 
autocorrelation function hypotheses comprises a plurality of samples and 
wherein the autocorrelation function for the estimate of the communications 
channel comprises a plurality of samples. 

10. A method according to Claim 1 wherein the communications 
channel comprises a radio channel. 

11. A method according to Claim 1 wherein the autocorrelation 
function comprises a plurality of samples and wherein generating the 
autocorrelation function comprises reducing a frequency error of the plurality 
of samples. 
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12. A method according to Claim 1 further comprising: 
providing a second estimate of the communications channel using the 

estimate of the Doppler spread for the communications channel. 

13. A method according to Claim 1 wherein providing an estimate of 
the communications channel comprises receiving a plurality of samples of 
data over the communications channel wherein the samples of data are used 
to generate the channel estimates, the method further comprising: 

5 performing a cyclical redundancy check on the samples of data used to 

generate the channel estimates; 

updating a long term Doppler estimate with the estimate of the Doppler 
spread when the samples of data pass the cyclical redundancy check. 

14. A method according to Claim 1 wherein providing an estimate of 
the communications channel comprises: 

receiving a plurality of pilot symbols over the communications channel; 
reducing phase ambiguity between the pilot symbols; and 
5 using the pilot symbols having reduced ambiguity to provide the 

estimate of the communications channel. 

15. A Doppler spread estimator for estimating a Doppler spread for 
a communications channel, the Doppler spread estimator comprising: 

a plurality of autocorrelation function hypotheses corresponding to a 
respective plurality of Doppler spread estimate hypotheses; 

a channel estimator that estimates the communications channel; 
an autocorrelation generator that generates an autocorrelation function 
for the estimate of the communications channel; 

an autocorrelation function hypothesis tester that selects one of the 
autocorrelation function hypotheses to approximate the autocorrelation 
function for the estimate of the communications channel; and 

a Doppler spread hypothesis selector that selects one of the Doppler 
spread estimate hypotheses corresponding to the selected autocorrelation 
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function hypotheses as an estimate of the Doppler spread for the 
communications channel. 



16. A Doppler spread estimator according to Claim 15 wherein the 
hypothesis tester compares the autocorrelation function for the estimate of the 
communications channel with each of the plurality of autocorrelation function 
hypotheses. 

17. A Doppler spread estimator according to Claim 15 wherein the 
hypothesis tester selects one of the plurality of autocorrelation function 
hypotheses most closely approximating the autocorrelation function for the 
estimate of the communications channel. 

18. A Doppler spread estimator according to Claim 15 wherein the 
hypothesis tester generates a plurality of error signals respectively 
corresponding to the plurality of autocorrelation function hypotheses wherein 
each of the error signals represents a difference between the respective 

5 autocorrelation function hypothesis and the autocorrelation function of the 
estimate of the communications channel, and wherein the hypothesis tester 
compares the error signals to select the autocorrelation hypothesis to 
approximate the autocorrelation function for the estimate of the 
communications channel. 

19. A Doppler spread estimator according to Claim 18 wherein the 
hypothesis tester chooses an error signal representing a least difference 
between the corresponding autocorrelation function hypothesis and the 
autocorrelation function for the estimate of the communications channel. 

20. A Doppler spread estimator according to Claim 19 wherein the 
hypothesis tester chooses a least one of the error signals. 

21 . A Doppler spread estimator according to Claim 1 8 wherein the 
hypothesis tester averages each of the error signals to provide averaged error 
signals wherein the hypothesis tester compares the averaged error signals. 
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22. A Doppier spread estimator according to Claim 21 wherein the 
hypothesis tester averages each of the error signals using one of block 
averaging, running averaging, and sliding window averaging. 

23. A Doppier spread estimator according to Claim 15 wherein each 
of the plurality of autocorrelation function hypotheses comprises a plurality of 
samples and wherein the autocorrelation function for the estimate of the 
communications channel comprises a plurality of samples. 

24. A Doppier spread estimator according to Claim 15 wherein the 
communications channel comprises a radio channel. 

25. A Doppier spread estimator according to Claim 15 wherein the 
autocorrelation function comprises a plurality of samples and wherein the 
autocorrelation generator reduces a frequency error of the plurality of 
samples. 

26. A Doppier spread estimator according to Claim 15 wherein the 
channel estimator provides a second estimate of the communications channel 
using the estimate of the Doppier spread for the communications channel. 

27. A Doppier spread estimator according to Claim 15 wherein the 
channel estimator receives a plurality of samples of data over the 
communications channel wherein the samples of data are used to estimate 
the communications channel, wherein the channel estimator performs a 

5 cyclical redundancy check on the samples of data used to generate the 
channel estimates, the Doppier spread estimator further comprising: 

a long term Doppier estimator wherein the long term Doppier estimator 
is updated with the estimate of the Doppier spread when the samples of data 
pass the cyclical redundancy check. 

28. A Doppier spread estimator according to Claim 15 wherein the 
channel estimator receives a plurality of pilot symbols over the 
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communications channel, reduces phase ambiguity between the pilot 
symbols, and uses the pilot symbols having reduced ambiguity to provide the 
5 estimate of the communications channel. 

29. A method for receiving communications, the method comprising: 
receiving a signal over a communications channel wherein the signal 

represents data from a remote transmitter; 

generating an estimate of the communications channel over which the 
5 signal is received; 

generating an autocorrelation function for the estimate of the 
communications channel; 

selecting one of a plurality of autocorrelation function hypotheses 
corresponding to a respective plurality of Doppler spread estimate hypotheses 
10 to approximate the autocorrelation function for the estimate of the 
communications channel; 

selecting one of the Doppler spread estimate hypotheses 
corresponding to the selected autocorrelation function hypotheses as an 
estimate of the Doppler spread for the communications channel; and 
15 reproducing an estimate of the data transmitted by the remote 

transmitter. 

30. A method according to Claim 29 wherein selecting one of the 
autoconnection function hypotheses comprises comparing the 
autocorrelation function for the estimate of the communications channel with 
each of the plurality of autocorrelation function hypotheses. 

31 . A method according to Claim 29 wherein selecting one of the 
autocorrelation function hypotheses comprises selecting one of the plurality of 
autocorrelation function hypotheses most closely approximating the 
autocorrelation function for the estimate of the communications channel. 

32. A method according to Claim 29 wherein selecting one of the 
autocorrelation function hypotheses comprises: 

21 



BNSDOCID: <WO 0113537A1_L> 



WO 01/13537 



PCT/USOG/21081 



generating a plurality of error signals respectively corresponding to the 
plurality of autocorrelation function hypotheses wherein each of the error 
5 signals represents a difference between the respective autocorrelation 
function hypothesis and the autocorrelation function of the estimate of the 
communications channel; and 

comparing the error signals to select the autocorrelation hypothesis to 
approximate the autocorrelation function for the estimate of the 
10 communications channel. 

33. A method according to Claim 32 wherein comparing the error 
signals comprises choosing an error signal representing a least difference 
between the corresponding autocorrelation function hypothesis and the 
autocorrelation function for the estimate of the communications channel. 

34. A method according to Claim 33 wherein choosing the error 
signal comprises choosing a least one of the error signals. 

35. A method according to Claim 32 wherein comparing the error 
signals is preceded by: 

averaging each of the error signals to provide averaged error signals 
wherein comparing the error signals comprises comparing the averaged error 
5 signals. 

36. A method according to Claim 35 wherein averaging each of the 
error signals comprises one of block averaging, running averaging, and sliding 
window averaging. 

37. A method according to Claim 29 wherein each of the plurality of 
autocorrelation function hypotheses comprises a plurality of samples and 
wherein the autocorrelation function for the estimate of the communications 
channel comprises a plurality of samples. 

38. A method according to Claim 29 wherein the communications 
channel comprises a radio channel. 
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39 A method according to Claim 29 wherein the autocorrelation 
function comprises a plurality of samples and wherein generating the 
autocorrelation function comprises reducing a frequency error of the plurality 
of samples. 

40. A method according to Claim 29 further comprising: 
providing a second estimate of the communications channel using the 

estimate of the Doppler spread for the communications channel. 

41 . A method according to Claim 29 wherein providing an estimate 
of the communications channel comprises receiving a plurality of samples of 
data over the communications channel wherein the samples of data are used 
to generate the channel estimates, the method further comprising: 

5 performing a cyclical redundancy check on the samples of data used to 

generate the channel estimates; 

updating a long term Doppler estimate with the estimate of the Doppler 
spread when the samples of data pass the cyclical redundancy check. 

42. A method according to Claim 29 wherein the data received from 
the remote transmitter comprises pilot symbols, and wherein generating an 
estimate of the communications channel comprises: 

reducing phase ambiguity between the pilot symbols; and 
5 using the pilot symbols having reduced ambiguity to provide the 

estimate of the communications channel. 

43. A receiver comprising: 

a radio receiver and converter that receives a signal over a 
communications channel; 

a channel estimator that estimates the communications channel 
5 responsive to the signal received over the communications channel; 

a plurality of autocorrelation function hypotheses corresponding to a 
respective plurality of Doppler spread estimate hypotheses; 
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an autocorrelation generator that generates an autocorrelation function 
for the estimate of the communications channel; 

an autocorrelation function hypothesis tester that selects one of the 
autocorrelation function hypotheses to approximate the autocorrelation 
function for the estimate of the communications channel; and 

a Doppler spread hypothesis selector that selects one of the Doppler 
spread estimate hypotheses corresponding to the selected autocorrelation 
function hypotheses as an estimate of the Doppler spread for the 
communications channel. 

44. A receiver according to Claim 43 wherein the hypothesis tester 
compares the autocorrelation function for the estimate of the communications 
channel with each of the plurality of autocorrelation function hypotheses. 

45. A receiver according to Claim 43 wherein the hypothesis tester 
selects one of the plurality of autocorrelation function hypotheses most closely 
approximating the autocorrelation function for the estimate of the 
communications channel. 

46. A receiver according to Claim 43 wherein the hypothesis tester 
generates a plurality of error signals respectively corresponding to the plurality 
of autocorrelation function hypotheses wherein each of the error signals 
represents a difference between the respective autocorrelation function 
hypothesis and the autocorrelation function of the estimate of the 
communications channel, and wherein the hypothesis tester compares the 
error signals to select the autocorrelation hypothesis to approximate the 
autocorrelation function for the estimate of the communications channel. 

47. A receiver according to Claim 46 wherein the hypothesis tester 
chooses an error signal representing a least difference between the 
corresponding autocorrelation function hypothesis and the autocorrelation 
function for the estimate of the communications channel. 
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48. A receiver according to Claim 47 wherein the hypothesis tester 
chooses a least one of the error signals. 



49. A receiver according to Claim 46 wherein the hypothesis tester 
averages each of the error signals to provide averaged error signals wherein 
the hypothesis tester compares the averaged error signals. 

50. A receiver according to Claim 49 wherein the hypothesis tester 
averages each of the error signals using one of block averaging, running 
averaging, and sliding window averaging. 

51 . A receiver according to Claim 43 wherein each of the plurality of 
autocorrelation function hypotheses comprises a plurality of samples and 
wherein the autocorrelation function for the estimate of the communications 
channel comprises a plurality of samples. 

52. A receiver according to Claim 43 further comprising: 

a signal processor coupled between the Doppler spread estimator and 
the channel estimator wherein the signal processor modifies«the operation of 
the channel estimator responsive to the estimate of the Doppler spread for the 
5 communications channel. 

53. A receiver according to Claim 43 wherein the communications 
channel comprises a radio channel. 

54. A receiver according to Claim 43 wherein the autocorrelation 
function for the estimate of the communications channel comprises a plurality 
of samples and wherein the autocorrelation generator reduces a frequency 
error of the plurality of samples. 

55. A receiver according to Claim 43 wherein the channel estimator 
provides a second estimate of the communications channel using the 
estimate of the Doppler spread for the communications channel. 
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56. A receiver according to Claim 43 wherein the signal received 
over the communications channel comprises a plurality of samples of data 
wherein the channel estimator uses the samples of data to estimate the 
communications channel, wherein the channel estimator performs a cyclical 

5 redundancy check on the samples of data used to generate the channel 
estimates, the receiver further comprising: 

a long term Doppler estimator wherein the long term Doppler estimator 
is updated with the estimate of the Doppler spread when the samples of data 
pass the cyclical redundancy check. 

57. A receiver according to Claim 43 wherein the signal received 
over the communications channel comprises a plurality of pilot symbols, 
wherein the channel estimator reduces phase ambiguity between the pilot 
symbols, and wherein the channel estimator uses the pilot symbols having 

5 reduced phase ambiguity to provide the estimate of the communications 
channel. 



_0113537A1_I_> 



26 



WO 01/13537 



1/4 



PCTAJS00/21081 



3 I 



TRANSMITTER 




FIG. 1. 




^23o 

RADIO 
RECEIVER 

AND 
CONVERTER 



25a 



1 



Ra 



CHANNEL 
ESTIMATOR 



SYMBOL 
ESTIMATOR 



29a 
V 




0113537A1_L> 



WO 01/13537 



2/4 



PCT/US00/21081 




TRANSMITTER 



FIG. 3. 



I 



Rb 



RADIO 
RECEIVER 
-ANB- 



CONVERTER 



SYMBOL 



CHANNEL 

ESTlMAreR-F=l-ESTIMATOR-|^- 



29b 



DOPPLER 
SPREAD 
ESTIMATOR 



SIGNAL 
PROCESSOR 

^31 




TRANSMITTER 



FIG. 4. 



RADIO 
RECEIVER 

AND 
CONVERTER 



I 



Rc 



25c 



CHANNEL 
ESTIMATOR 



27c^ 



T 



DOPPLER 
SPREAD 
ESTIMATOR 



SIGNAL 
PROCESSOR 



) 0113537A1_L> 



WO 01/13537 



PCT/USOO/21081 



3/4 



CHANNEL 
ESTIMATES 



DOPPLER SPREAD 
ESTIMATOR 



FIG. 5. 



AUTO CORRELATION 
CALCULATOR 



51 



r 



53 















Hi 




H 2 




H 


3 


• • • 


"K 



ej (m) $2 l™) e 3 ( m ) 



e K (m) 



AVERAGER 



'av,l e av,2 e av,3 



e av,K 



MINIMUM ERROR SELECTOR 



57 



DOPPLER SPREAD HYPOTHESIS SELECTOR 



59 



DOPPLER 



VO 0113537A1_I_> 



WO 01/13537 



4/4 



PCT/US00/21081 



1 


„_^_^ Hj@20Km/hr 


^ 


\ ^v^^\H2@40Km/hr 


C3 


\ He® \ \ 
\100Whr ^0Krn/hr\H 3 @60Km/br" 


0 






— ^*<^S^^^— t 




FIG. 6. 



H, @20KmAr 


R H ,i(0) 






Rh,i(t z ) 


H 2 @40Km/hr 


R Hf2 (0) 


R H ,2(T,) 




" H,2<T Z ) 


H 3 @60Km/hr 


R H>3 (0) 


R H/3 (Ti) 




R H f 3<V 


H 4 @80Ktn/hr 


R H(4 (0) 


R H,4^i) 




R H,4^ z ) 


H 5 @1000Km/hr 


R H(5 {0) 






R H,5<V 



FIG. 7. 



_0113537A1J_> 



INTERNATIONAL SEARCH REPORT 



Intern -»nal Application No 

PCT/US 00/21081 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 H04B7/01 



According to International Patent Classification (IPC) or 10 both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system toil owed by classification symbols) 

IPC 7 H04B 



Documentation searched other than minimum oocumentaiion to the extent that such documents are included in the fields searched 



Electronic data base consulted during the wterruumai search (name of data base and, where practical search terms used) 



PAJ, EPO-Internal , WPI Data, INSPEC 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category • Citation of document, wdh w^on. 



Appropriate, of the relevant passages 



Relevant to claim No. 



p,x 



DATABASE WPI 

Section EI, Week 199945 

Derwent Publications Ltd., London, GB; 

Class W02, AN 1999-535416 

XP002152689 

& JP 11 234190 A (OKI ELECTRIC IND CO LTD) 

, 27 August 1999 (1999-08-27) 

abstract 



1,15,29, 
43 



US 5 606 580 A (MOUROT CHRIST0PHE 
25 February 1997 (1997-02-25) 
abstract 

claim 1 

figure 4 



ET AL) 



1,15,29, 
43 



-/— 



Further documents are listed in the 



of bos C 



ID 



Paten) famHy members are listed in annex. 



■ Special categories of cited documents : 

"A" document denning the general siaie oi in# *rt ***** a not 

considered to be of particular retevanc* 
'E' earlier document but published on or after ttm mttwvi *nnai 

filing date 

"L" document which may throw doubts on pnonry ci^to m or 
which is cited to establish the publication o#e ct mooter 
citation or other special reason (as spooled) 

•C document referring to an oral disclosure, use. e*tt*x*on or 
qther means 

'P' document published prior to the international hfcng date but 
later than the priority date claimed 



later document published after the international filing date 
or pnority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

document of particular relevance: the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

document of particular relevance; the claimed invention 
cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person slatted 
in the an. 

document member of the same patent family 



Date of the actual completion of the international saarcn 



10 January 2001 



Name and mailing address of the ISA 

European Pateni Office. P.8. 58ia Patemtaan 2 
NL - 2280 HV Rijswijk 
TeL (+31-70) 340-2040. Tx. 31 651 epo nl. 
Fax (+31-70) 340-30 16 



Date of mailing of the international search report 



18/01/2001 



Authorized officer 



Lustrini , D 



Forni PCT/1SA/210 (second sheet) (July 1992) 



page 1 of 2 



BNSDOCID: <WO 0113537A1_I_> 



INTERNATIONAL SEARCH REPORT 



Interr. >nal Application No 

PCT/US 00/21081 



C.(ContlnuatJon) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * Citation of document, with indication. where appropriate, of the relevant passages 



Relevant to claim No. 



EP 0 363 640 A (INT STANDARD ELECTRIC 
CORP) 18 April 1990 (1990-04-18) 
abstract 

page 8, line 25 - line 33 
page 10, line 18 - line 45 
figure 9 



1,15,29, 
43 



Form PCT/1SA/210 (continuation of second sheet) (July 1992) 



page 2 of 2 



BNSDOCID: <WO 0113537A1_1_> 



INTERNATIONAL SEARCH REPORT 

...formation on patent family member* 


Inter?. -naJ Application No 

PCT/US 00/21081 


Patent document 
cited in search report 


Publication 
date 


Patent family 
member(s) 


Publication 
date 



JP 11234190 A 27-08-1999 NONE 



US 5606580 A 25-02-1997 FR 2708162 A 27-01-1995 

EP 0635949 A 25-01-1995 

JP 7154361 A 16-06-1995 

NO 942704 A 23-01-1995 



EP 0363640 A 18-04-1990 US 4914699 A 03-04-1990 

DE 68927658 0 27-02-1997 
DE 68927658 T 24-07-1997 



Form PCT/1SA/210 (patem family amtaL) (July 1B92) 
BNSDOCID: <WO 0113S37A1 I > 



